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Abstract. The statistical properties of non-interacting bosons and fermions confined 
in trapping potentials are most easily obtained when the system may exchange 
energy and particles with a large reservoir (grand-canonical ensemble). There are 
circumstances, however, where the system under consideration may be considered as 
being isolated (micro-canonical ensemble). This paper first reviews results relating to 
micro-canonical ensembles. Some of them were obtained a long time ago, particularly 
by Khinchin in 1950. Others were obtained only recently, often motivated by 
experimental results relating to atomic confinement. A number of formulas are 
reported for the first time in the present paper. Formulas applicable to the case 
where the system may exchange energy but not particles with a reservoir (canonical 
ensemble) are derived from the micro-canonical ensemble expressions. The differences 
between the three ensembles tend to vanish in the so-called Thermodynamics limit, 
that is, when the number of particles and the volume go to infinity while the particle 
number density remains constant. But we are mostly interested in systems of moderate 
size, often referred to as being mesoscopic, where the grand-canonical formalism is not 
applicable. The mathematical results rest primarily on the enumeration of partitions 
of numbers. 
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1. Introduction 

Let us first recall that in Quantum Statistical Mechanics isolated systems should be 
treated in the microcanonical formalism, systems that may exchange energy with a 
heat bath in the canonical formalism and systems that may exchange both energy and 
particles with a bath in the grand-canonical formalism. In the latter case simple formulas 
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are easily obtained. For large systems the three ensembles are generally eqnivalent so 
that the grand-canonical formalism suffices for most applications. However, there is 
presently much interest in systems consisting of only a few particles, bosons or fermions. 
Recent experimental observations of atomic Bose-Einstein condensation motivated much 
theoretical work relating to mesoscopic devices. The fluctuation of the number of bosons 
occupying the ground level of harmonic traps was particularly considered. In such cases 
the grand-canonical ensemble is inappropriate and one must consider instead the less- 
easily handled micro-canonical and canonical formalisms. 

Many results about canonical and microcanonical ensembles are known but they 
are scattered in the literature and often require intricate mathematical methods. The 
present paper provides an essentially exhaustive and mathematically simple derivation 
of the statistics of non-interacting particles in the three ensembles. Fermions and bosons 
are treated on the same footing but particles such as photons and phonons that do not 
carry generalized charges are considered in a separate section. In our approach system 
energy levels and their degeneracies are considered. A short study of compound systems 
helps clarify how degeneracy occurs. The mathematical method rests on exact results 
relating to the microcanonical statistics, canonical and grand-canonical expressions 
being subsequently derived through averaging. This line of attack is often evoked as 
being the most natural one but it has apparently not been pursued in detail before. 
Microcanonical statistics was treated by Khinchin in 1950 but other statistics and 
degeneracies were not considered at that time. The microcanonical theory rests on the 
assumptions that the particles are undistinguishable and interact only weakly so that 
one-particle energies add up. The number of microcanonical conhgurations is assumed 
to be hnite. 

The purpose of this paper is not purely educational. Although the mathematical 
formalism requires only elementary combinatorics and the handling of generating series, 
the method naturally gives rise to several new results. This is the case in particular 
for formulas that facilitate estimations of certain quantities relevant to the canonical 
and grand-canonical ensembles. Complete results are given in the important case 
of evenly-spaced non-degenerate levels, one-dimensional harmonic potentials being a 
basic example. Recursive relations for occupancies are obtained that generalize results 
recently reported by Schonhammer |Q. It has been noticed long ago that microstates are 
in that case described by the partition of integers P, It is also known that hard-core 
interacting bosons in harmonic potentials may be treated as non-interacting fermions 
Furthermore, the micro-canonical statistics of a band of states is particularly relevant 
to spin systems P and semiconductor devices. 

The paper is essentially self-contained. Various formulas about occupancies, 
statistical weights, and partition functions for bosons and fermions are obtained in 
Sections 2 to 4, beginning with the microcanonical and ending up with the grand- 
canonical ensemble. Chargeless particles are treated in Section 5, and Section 6 
addresses issues relating to compound systems and degeneracy. Section 7 deals with 
the case of evenly-spaced levels. Relations between fermionic and bosonic statistics are 
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given in Appendix A, links between von Neumann-Shannon entropy and the classical 
formulations of entropy in Thermodynamics in Appendix B. 

2. Microcanonical statistics 

Consider an isolated system comprising N undistinguishable particles. The energy Sg of 
a particle in (quantum) state s is assumed to be non-negative without loss of generality. 
In what follows S denotes the set of states and E the set of energy levels. The degeneracy 
Qg of an energy level £ is the number of states corresponding to that level. For example 
the energy of a state s = (fci, ..., kd), where ki,... ,kd are non-negative integers, of a 
d-dimensional isotropic harmonic oscillator is given by the formula £ = hulk d/2) 
where hcu is a constant and k = ki kd- The degeneracy of the energy level e 

is equal to the number of ways of obtaining k by summing up d nonnegative integers, 
namely Qe = {k + d — 1)\/{k\{d — 1)!). 

A system is characterized by N and the total energy U it contains. Because of our 
assumption that the one-particle energies add up, U is the sum of N £-values (e G E), 
some of them possibly occuring more than once. In the sequel U denotes the set of all 
hnite sums of elements of E, with the convention that U contains 0. Given N E N and 
t/ G U several conhgurations may occur when U is a sum of N values in E. Letting Ug 
denote the number of particles in state s, a (A^, 17)-microstate is a family u = {i'g)g£s 
of non-negative integers (in other words, z/ is a mapping from S to M) fullhlling the 
conditions 

'^Ug = N, and ^ UgSg = U. (1) 

sgS sgS 

As an example consider a system of particles whose energy levels are £ = 2A; — 1 -|-/^, 
with k, I positive integers (this is the case for charged single-spin particles in uniform 
magnetic helds under some conditions). Figure |I] shows the 22 partitions of t/ = 22 
among N = 3 particles. Because of the degeneracies some of the partitions shown 
are associated with several microstates. The number Vg of bosons in some state s is 
unconstrained while the numbers of fermions in a state may only be 0 or 1. Accordingly, 
the hrst marked (|) partition corresponds to two bosonic microstates since the £=18 
level is twice degenerate and to zero fermionic microstate since the e=2 level is not 
degenerate. The second marked partition in Fig. |l| may correspond to a fermionic 
microstate because two distinct states have energy £=10, namely (i.e., k=l, 1=3) 
and S5,i (i.e., k=5, 1=1), or to three bosonic microstates where the two-particles states 
for £=10 are respectively: (51,3,55,1), (51,3,51,3) and (55,1,55,1). 

The system statistical weight W{N,U), i.e., the number of distinct {N,U)- 
microstates, is assumed to be hnite. This implies hnite degeneracies. Indeed, W{1,U) 
is zero if U is not a one-particle energy level, otherwise it is the degeneracy at this level. 
Notice that since 0 G U by convention hF(0, 0) = 1. Summing up the numbers on top of 
Fig. g we see that W{/3, 22) is equal to 34 for bosons and 21 for fermions. The purpose 
of formulas given below is to provide such numbers directly. 
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Figure 1. This figure represents the (TV = 3, U = 22)-microstates: The energy levels 
e and their degeneracy are shown on the left. At the top of the figure the number 
of microstates corresponding to each partition is given for bosons (b) and fermions (f). 


The microcanonical occupancy N^{N, U) of energy level £ is defined as the average 
number of particles at this energy level, the W{N, U) microstates being equally likely. 
(If W{N, (7) = 0 we have N^iN, U) = 0). In Fig. |l| the occupancy of level £ = 10 is 
12/34 for bosons and 6/21 for fermions. Letting f2(iV, U) denote the set of all (iV, U)- 
microstates the occupancy may be written 

= p) 

^ ^ S,es=e ^£n(N,U) 

the hrst summation index s running through the states of energy £. Conditions (|l]) and 
dehnition (||) result in 

Ns{N, U) = N, and ^ N,{N, U)e = U. (3) 

£GE eSE 


For example, the hrst relation in (|^) derives from 

w{N.u)Y,r^.{N,u)^J2Y. E E "•= E E^-' 

eGE eGE s,£s=e uen{N,U) sGS iy£n{N,U) uGn{N,U) sGS 


Notice that at most NW{N, U) summands contribute to each of the above sums. In 
the sequel the symbol + refers to bosons and the symbol — to fermions. The symbol ± 
should be replaced by + for bosons and by — for fermions. 

The present article rests mainly on the following recursion relation, to be proven 
later in this section 


N^{N,U) 


W^{N -1,U -e) 
hF±(Af, U) 


(^7,±iV±(iV 


l,U-e)). 


( 4 ) 
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Since W{N — n,U — ne) vanishes for n > N, iterating N times the above formula yields 

N 

Nf(N, U) = -n,U- ns). (5) 

Combining and provides recursive formulas for the statistical weight 

1 ^ 

W^{N, U) = — EE(±ir ^QeW^iN -n,U -ne) (6) 

n=l esE 
1 ^ 

= ^EE(±ir ^egeW^{N -n,U -ne). (7) 

n=l eSE 

For example, the hrst one obtains as follows: 

N 

NW^{N, U) = U)W^{N, U) = ^ - n, t/- ne). 

eGE eGE n=l 

Notice that the sums above involve only the values of e for which N^{N, U) ^ 0, that is, 
at most NW{N, U) e-values. Formula (|^) appears in |]^, § in a different form (summation 
is done on one-particle states instead of system energy levels) for bosons, and in 0 for 
fermions without degeneracy. Although less interesting than for a direct computation 
of W{N, U), formula (|^) will prove useful later. 

To establish (|^) consider a state s. For bosons m shows that the relation Vg = u'^ — l 
dehnes a one-to-one correspondence between the set of [N — 1, 7/— £)-microstates u and 
the set of {N, t/)-microstates z/' with > 0. Thus 

5^ (l + i^.) = WnN-l,U-e)+ Y. (8) 

y'eei+(N,U) um+{N-l,U-e) ueQ.+ {N-l,U-e) 

The above result amounts to removing one particle in state s from each (A^, U)- 
microstate that exhibits such a particle and counting the remaining particles in state s. 
For fermions the number of {N, 17)-microstates with one particle in state s equals the 
number of [N — 1, f/ — £)-microstates with no particle in state s, that is, all of them 
except those having one particle in that state. Thus 

5^ C, = W-(N-l,U-e)- Y (9) 

u'eei-(N,u) uen-{N-i,u-£) 

By (H), summing extremal members in (j|) and (§) for all states s of energy e the result 
in (^) is obtained. 

Formulas (|^) and have apparently not been published before in their general 
form. Related expressions are given by Khinchin |l| under the assumption of evenly- 
spaced non-degenerate energy levels. Khinchin accounts for degeneracy by repeating 
energy levels. This amounts to consider quantum states instead of energy levels 0]. For 
bosons Eq. (|^) has appeared in a different form in [^. For the case of non-degenerate 
fermionic systems the formulas were given in [^]. 

It is thus clear that microcanonical statistics rely on the sole knowledge of statistical 
weights, the computation of which requiring only the energy-levels values and their 
degeneracy. 
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When the system may exchange energy with a medinm of arbitrarily large size at 
temperatnre T (heat bath), the probability that a particnlar (A^,M)-microstate occnrs 
at eqnilibrinm is proportional to where 


q = exp(- 


k^T 


), 


according to Boltzmann and Gibbs. The partition function 


Z{N,q) = J2w{N,u)q^ 

u£V 


( 10 ) 


normalizes the probability. Notice that, even thongh all snms are finite in the 
previous section, the summation index u takes on infinitely many contributing values. 
Summability requires that the energy levels £ form a finite or unbounded countable set. 
If the levels are discrete this means that they must form a finite or unbounded increasing 
sequence. Convergence of the partition function occurs if, and only if, Z(l, q) converges. 
Indeed, 0 < Z~{N,q) < Z+(iV, g) < Z{l,q)^. We assume in the sequel that Z{l,q) 
converges. 

The system is characterized by the number N of particles and the temperature T, 
or N and q equivalently. The average energy U{N, q) follows from 

U{N, q) = Y. ^)- (11) 


The system may be characterized by N and U since U{N, q) is an increasing function 
of q. Indeed, the energy variance, a positive quantity, is given by 

Var(G) = —^ u{u - 1)W{N, u)q^ + f/(iV, q) - U{N, qf 

^ ueu 

«>■ 

Energy variance is therefore linked to heat capacity c{N, T) = dU{N, q)/dT by Var(f/) = 
k-QT‘^c{N,T). It is worthwhile to point out that two systems in the canonical ensemble 
have the same heat capacity if, and only if, their partition functions coincide except 
for a multiplicative factor of the form ag^, where a and b are real numbers that do not 
depend on T. This result follows easily from (|^. 

The canonical occupancy N^^N^q) at level £ is the average number of particles at 
this level. It may therefore be calculated by averaging the microcanonical occupancies 
at level e with W{N,u)q^ as a weight. Accordingly, relations result in 

q) = iV, and q)e = U{N, g). 

eeE eeE 


( 12 ) 
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Since W{N, u) vanishes for m < 0, it is not difficult to see ||] with the help of formula 
that 


N^{N,q) = q^ 


Z±(iV,g) 


{g,±Nf{N-l,q)), 


9s 


N 


ZHN,q) 


5^(±ir-'g-Z±(iV-n,g). 


(13) 


(14) 


n=l 


Equation (|1^) is derived by iterating (|T^) or averaging (^. A relation similar to (^) 
is given in |]in for the case of non-degenerate energy levels. But the assumption made 


in that paper that Nf^{N, q) and Nf^{N — l,q) are identical is valid for large A^-values 
only. The relations ([T^ ) and (^ were given in [pT[| . 

In general, a closed form expression of Z{N, q) has not been established, but this 
quantity can be evaluated recursively. Using 1U(0,0) = 1 and 1U(0, m) = 0 if m > 0, we 
obtain from (|6|) 


Z(0,g) = l, and Z(l,q} = 


(15) 


eSE 


Summing up (|T^) for £ running through E and making use of the hrst formula in ([I^) , 
we obtain for > 1 

a N 


Z*( V = ^(±l)"-‘Z±(iV - n, q)Z(l, q' 


(16) 


n=l 


This formula may also be derived directly from (|^). Alternative proofs are given in 
. Making use of (^ instead provides a similar formula for the computation 


of Z'(N,q) = ■^Z{N,q). Comparing this formula with the derivative of (pffi[) results in 
an other expression. These new formulas are respectively 


N 


Z'^{N,q) = J2(±qr-'Z^{N - n,q)Z'{l,q’‘) 


(17) 


n=l 


1 

N- 1 


N-l 

J2i±^r-^Z'^{N -n,q)Z{l,q-). 

n=l 


(18) 


The above relations (^) and (^) allow, in turn, the numerical computation of U{N, q) 
with the help of ([TlD. 

Note incidentally that Z{l,q) depends only on the energy levels and their 
degeneracy, not on the nature of the particles. Moreover, formulas ( 0 ) and ( 0 ) 
show that the main difficulty in the computation of Z{N,q) and N^{N,q) rests on 
the numerical evaluation of Z{l,q"'), that is, of the one-particle partition function at 
temperature T/n, for n going from 1 to N. There is a formula giving Z{N, q) in terms 
of Z{l,q) only: 


Z±(iV,g) 


E n 

Ufe)fc>i6P(Ar)fc>l 


(±l)^-iZ(l,g^)^fc 

z/fc! 


(19) 


f Index shifts in the calculation, such as X]ugu^(-^ — 1,u — e)g" = — l,u)q^, are 

justified by the following remark. Let u' G U. If there exists N d'H* such that W{N^ u — u') 0, then 

u G U if, and only if, u — u' G U. 
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where P(A^) stands for the set of unrestricted partitions of N, corresponding to the 
non-degenerate bosonic {N, f/)-microstates with E = N and U = N] see, e.g., 0. This 
formula, which may be established by induction on N (the proof is tedious but obvious), 
is a special case of more general relations that are not relevant to this paper. Related 
formulas are given in |^, and linked to symmetric polynomials in . 


4. Grand-canonical statistics 

If the system may exchange energy and particles with a bath at temperature T and 
fugacity z, the probability that a (n, n)-microstate occurs is proportional to The 

grand-partition function 

9 ) = u)z^q^ 

n>0 uGU 

that normalizes the probability reads 

1 


2+(z.9) = n7TV 


(1 — zq^Y^ ’ 




£GE ' eGE 

These formulas can be recovered by averaging both sides of 
and using (|^ so that 

n 

nZ^{n, q)z'^ = Y “ b ?) 3e(f' 

n>0 i=l eGE 


zqy\ (20) 

with z^ as a weight 


n>0 


£:GE n >0 i=l 

- E sdE(±‘«')“) (E i)r 


eGE 


n>l 


= z 




n>0 

Tq" 


1 =F zq’^ 


n>0 sEE 

Dividing both sides by 2 : X]n>o recognize that the two sides of each equality 

shown in O have the same logarithmic derivative with respect to 2 ;, so that they are 
proportional to one another. It follows that they are identical since they have the same 
constant term 1. 

The grand-canonical occupancy at level £ is calculated by averaging the 
microcanonical occupancy N^{n,u) with z‘^q'^W{n,u) as a weight, that is, by averaging 
the canonical occupancy N^n, q) with Z(n, q)z^ as a weight. Using ([I^ and introducing 
the chemical potential /i = k^T In ; 2 , we obtain the celebrated Bose-Einstein and Fermi- 
Dirac distributions 


Nf-{z, q) = 


9e 


Ne (ug) = — 


9e 


( 21 ) 


ij'-'-i •' r-‘ + i 

It follows from lnZ^(z, g) = TX]eGEg£^^(^ T zqY that the average number of 
particles and average energy are respectively 


d 




rAT 


— 


d 
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The expansion of ln(l =|= zq^) in power series of provides with an expression of 
In Z{z,q) in terms of the one-particle canonical partition fnnction, namely 

n 

In Z±(^, g) = V(±1)”->Z(1, (23) 

n 

n>l 

Accordingly, the average nnmber of particles and energy may be compnted as 

N{z,q) = J2i±^r~'Z{l,q-)z-, ( 24 ) 

n>l 

(25) 

n>l 


5. Photons and phonons. 

Particles that do not bear generalized charge snch as photons or phonons can be created 
or annihilated withont violating conservation rnles. Accordingly, the nnmber N of 
snch particles is not conserved in the micro-canonical and canonical ensembles. A 
f/-microstate is now a family u = {i's)s&§ of non-negative integers satisfying the sole 
condition 

= U. ( 26 ) 

sgS 

The micro canonical statistical weight W {U) is obtained by snmming np W (n, U) over 
all n values. Let eo = niinE be the lowest one-particle energy level. If £o > 0 then usq 
exceeds U at sufficiently large values of n so that W{U) is hnite. This arguments also 
applies to fermionic systems with eg = 0 because the number of fermions that occupy 
this level cannot exceed its degeneracy. But if Eq were equal to zero in bosonic systems 
the addition of a particle at level Eq would entail that W{n,U) < W{n -\- 1,U) and 
W{U) would be equal to zero or inhnity. We therefore assume that the ground-state 
energy of chargeless bosons is non-zero. No generality is lost since such particles cannot 
be detected at zero energy. 

The microcanonical occupancy of the energy level e is the average number of 
particles at level e. Arguments similar to the ones given in Section ^show that (||)-(||) 
may be written as 

Y,Ne{U)E = U, 

£GE 

NtiU) = ± NfiU - e)) ( 27 ) 

= Fl7m P8) 

' ' n=l 

Equation (|^) becomes 

. [U/s] 

E (±i)”-‘ir*(G -n£). 

£<U n=l 


( 29 ) 
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Canonical ensemble and grand-canonical ensemble with fugacity z = 1 are presently 
equivalent. Indeed, averaging as in §3 the above relations with as a weight gives 

Z(g) = 5^iy(w)g“, U{q) = q^\nZiq), iV,(g)£ = 17(g), (30) 




and 


= - 


9. 


- 1 


NrM = - 


eSE 


9e 


q-^ + l 


(31) 


Inserting the latter expressions in the third equation of ([30|) and integrating, the 
expression of InZ^(g) is obtained from the second equation of (|30[) yielding 

= Z-(9) = 11(1+«')’•■ (32) 


eGE ' ^ eSE 

If we expand ln(l ib g^) in power series of g^ as was done earlier in §4 we obtain 


ltiZ±(9) = 5^bt!£72(i,," 

c J T) 


(33) 


n>l 


The average energy may be derived from this expression according to (^) and the 
average number of particles is obtained from a summation of (^I]) and expansion of 
g^/(l ± q^). These two quantities are given by 

Af(9) = 5](±l)"-‘Z(1,9”), (7(9)=5](±1)”-‘Z'(1,9”)9” (34) 


n>l 


n>l 


6. Compound systems 

A compound system is a mixture of several sub-systems comprising undistinguishable 
particles, the particles of any two sub-systems being distinguishable. We assume that 
mixing does not affect the one-particle energy-level values and consider two sub-systems 
only, generalization to a greatest number of sub-systems being straightforward. 

Consider a system composed of two sub-systems with a total of W + -^2 = -^ 
particles and energy Ui + U 2 = U. If neither energy nor particles are allowed to 
be exchanged the number W{N, U) of conhgurations of the compound system is the 
product of the statistical weights of the sub-systems and the occupancies add up, that 
is 

W{N, U) = Wi(W, 17i)W2(iV2, t/ 2 ), (35) 

7V,(7V, U) = W,e(A^i, U,) + N2,e{N2, U^). (36) 

Suppose now that the number of particles in each sub-system is constant but that 
some energy may be transferred from one sub-system to the other. For example the two 
sub-systems may comprise electrons of spin 1/2 and -1/2 respectively, spin flip being not 
allowed. Each possible energy splitting between the two sub-systems has to be taken in 
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account. Since all the possible configurations are equally likely to occur any (iVi,Mi)- 
microstate of the first sub-system has a statistical weight W 2 {N 2 ,U — Ui). Therefore, 
the number of possible configurations of the system and the occupancy of level e are 


W{N, U) 
Ns{N, U) 


Y, bTl(iVl,Wl)1^2(iV2,W2), 

ui+U2=U 


ui-\-U 2 =U 


Wi{Nr,ui)W2{N2,U2) 

W{N, U) 


+ N2,e{N2,U2)). 


(37) 

(38) 


As one may expect for sub-systems that only exchange energy the canonical-ensemble 
formulas are simple. Remembering that at equilibrium the temperatures and chemical 
potentials of the two sub-systems equalize, the above expressions give 


Z{N,q) =Z,{N,,q)Z 2 {N 2 ,q), (39) 

N,{N, q) = iVi,,(iVi, q) + iV 2 ,,(iV 2 , q). (40) 


If the sub-systems may exchange particles, for instance if spin flip is allowed 
in a system of electrons, similar arguments show that the statistical weight and the 
occupancy of level e are 

W{N,U) = E E Wi{ni,ui)W2{n2,U2), (41) 

ni-\-n2=N ui-\-U2=U 


N,(N,U) 


E E 

ni-\-n2=N ui-\-U2=U 


Wi(ni,Ui)W2(n2,U2) 
W{N, U) 


Ml) -F N2,e{,n2, U2)). 


(42) 


Partition and grand-partition functions are easily derived: 
Z{N,q)= Y Zi{ni,q)Z 2 {n 2 ,q), 

ni-\-n2=N 

Z{z,q) = Zi{z,q)Z 2 {z,q). 


(43) 

(44) 


Such a compound systems is equivalent to a single system with the following 
characteristics: E is the union of the sets of energy levels of the subsystems and the 
degeneracies are added. Because of the possible exchange of energy and particles the 
two sub-systems become undistinguishable from one another. This fact is ascertained by 
the expression obtained for the grand-partition function of the compound system. But 
this conclusion does not apply when the sub-systems only exchange energy as illustrated 
in 0. 


7. Non-degenerate system with evenly-spaced levels 
7.1. Microcanonical relations 

When the energy levels are evenly spaced we may assume without loss of generality that 
the energy levels are non-negative integers. We may indeed take the energy-level spacing 
as the energy unit. On the other hand shifting the energy levels by one is tantamount 
to shifting the total energy of the system by the number of particles. The number of 
configurations and the occupancies are unaffected. 
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We first restrict our attention to a band of evenly-spaced one-particle energy levels, 
namely E = {0, 1,... ,B}. Let us recall that a partition of [/ G N is a non-increasing 
sequence of nonnegative integers summing up to U. If bosons are considered, a {N, U)- 
microstate is a partition of U with at most N parts none of them exceeding B. Thus 
the statistical weight is the number of such partitions. An example is given in Figure 


level 


6 






• 

• 






5 


• 






• 

• 




4 



• 






• • 

• 

• 


3 




• • 




• 


• 


• • 

2 




• 

• • 

• 



• 


• • 

• 

1 



• 

• • 

• 


• • 


• 

• 



0 

• •• 

• • 

• 

• 


• 


• 

• 





o 

II 



C/ = 5 





II 

00 





Figure 2. Bosons, {N = 3, {7)-microstates for B = 6. 


With a usual notation, the statistical weight 

W+{N,U) =p{B,N,U) (45) 

vanishes when U < 0 and U > NB. Several classical properties of the p{B, N,U)^s 
are given in m- For fermions, only the partitions of U with distinct parts are to be 
considered. Considering upward shifts from the smallest energy t4nin = N{N — l)/2 
and letting P = B — N 1 denote the maximum number of adjacent empty levels, it is 
easy to see that these partitions are the partitions of the added energy R = U — f/min 
into at most N parts none of them exceeding P, that is 


W-{N,U) =p{B - N + 1,N, U - N{N -l)/2) =p{P,N,R). (46) 

Figure ^ illustrates the above argument with reference to Figure The numbers in 
parentheses count the number of energy jumps from the initial positions and relate the 
fermionic (3,8)-microstates to the corresponding bosonic (3,5)-microstates. 

In the present situation of evenly-spaced levels further micro-canonical relations 
may be established. The following reasoning may be formalized by using the definition 
in (|l]) of microstates. Consider first a bosonic system consisting of N particles and 
energy U. Recall that W^{N — 1,U — e) is the number of {N, f/)-microstates with at 
least one particle at level e. There are therefore W~^{N, U) — W^{N — 1,U) microstates 
with no particle at level 0. Shifting downward the energy levels by 1 these microstates 
provide all the {N, U — Af)-microstates except those that have one particle at least at 
level B. In other words 


W+{N, U)-W+{N-1,U) = W+{N,U-N)-W+{N-1,U-N-B).{A7) 

Let us denote for short Mk(N,U) = W{N,U)Nk{N,U). Substituting N — n and 
U — n{k + 1) in place of N and U in (^^ and summing both sides with n running 
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Figure 3. Fermions, {N = 3,17)-microstates for B = 6. 


from 0 to N, we obtain from (|^) a relation involving the same nnmber N of particles 
but successive levels k,k 1. This relation reads for 0 < k < B 

W+(N,U) + M+^^(N,U)-M*^^{N,U + k+l) = 

W*(N,U - N)+M^(N,U - N)-Mt(N,U - N - B+ k). ' ’ 

Indeed, using (|^) we have successively 

W^{N -n,U -n{k + 1)) = iy+(iV, U) + iy+(iV, t/)iV++i(Ar, U), 

n>0 

Y -l-n,U -n{k + l)) = Y - n,U + k + I - n{k + 1 )) 

n>0 n>l 

= W+{N, U + k + l)iV++i(iV, U + k + 1), 

the right-hand terms in (^) being handled analogously. 

Consider now a fermionic system consisting of N particles and with total energy 
U. We have 

W-{N, U) + W-{N -1,U - B -1) = W-{N, U -N) + W-{N -1,U - N +1). (49) 

Indeed, the left-hand side of this expression counts the union of the (iV, 17)-microstates 
and the hctitious {N, 17)-configurations with 1 particle at the hctitious additional level 
B 1. Depending on whether the above {N, f/)-con£gurations have 0 particle at level 
0 or not, shifting downward the levels by 1 results either in a {N, U — iV)-microstate 
or in a {N,U — A^)-con£guration with 1 particle at the hctitious level —1. This gives 
the right-hand side. Notice further that analogously the occupancy of the lowest level 
satishes the relation 

Mo{N, U) + Mq{N -l,U -B-l) = W-{N -l,U - N + 1), (50) 

because shifting the levels yields all the [N, [/)-configurations with 1 particle at level 
— 1. In the same way as in (^Sf) one derives from (^), we obtain from (^^ for 0 < /c < 5 

W-[N,U)- A4-^,(V,(7) + (N,U-B + k) = 

W-{N,U -N) - h’q{N,U - N) + M;{N,U - N + l + k). 


( 51 ) 
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Replacing Mq{N, U) in (^OD by its value given by and substituting N to N — 1 we 
obtain 

U) = Mo{N, U-B-l) + W-{N, U) - W-{N, U - N). 

This relation shows that formula (|^) still holds for fc + 1 = 0. 


1.2. Canonical relations 


It is well known that the generating series for U of the numbers p{B, N, U) in 
the gaussian polynomial 

NB 

J2p{B,N,U)q^ 


IS 


(?) 


B+N 


U=0 


(?)b(?)v’ 


where {q)n is the commonly employed notation for the product (1 —g)(l — 
and (g)o = 1. Accordingly, 

{q)b+N ry-fT^T \ N(N-1)I2 (?)-B+l 


Z+(7V,g) = 


Z-{N,q) = 


( 1 -?") 


(52) 


(?)b(?)v’ ’ (?)B+l-v(?)Ar 

In order to keep the paper self-contained let us point out that the latter expressions 
may be easily derived from the previous microcanonical relations. According to (^) 
and (|4^) it suffices to consider Z~^{N,q). Firstly 


Z{Cq) = l + q + --- + q = 


1 




B+l 


Q 


Secondly, averaging (EZD with q^ as a weight and collecting the results we obtain the 
recurrence relation 

X — qB+n 

Z*(N,q)= Z*(N-l,q), 

from which the derivation of Z^{N.,q) in (^) is straightforward. 

Averaging (|48|) and (^) with q^ as a weight, dividing by Z{N, q) and multiplying by 
qk+i^ provides simple recursive formulas for the computation of canonical occupancies. 
These formulas are new for bosons, and generalize |@] for fermions. For t] < k < B they 
are 

(l-y+')Af++i(Af,,) =y+‘(l-«y-9"(y+'-9®+‘)A'^(A'.9), (53) 

(/+1 _ g»+>)iV-^,(iV,g) = y+‘(l - 9 ") - "("(I - (V 9)- (54) 

As noticed after (|^) the latter formula holds for A; -|- 1 = 0 so that 


No(N,q) = 


q 


N 




There is apparently no such simple expression for Nq {N, q). This quantity may be 
calculated from (1^) with the initial value Nq (1, g) = 1/Z(1, g). The expression obtained 
directly from (|T4|) is 

WN+B ^ l?jn 
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If we let -B —oo in the previous formulas we obtain expressions applicable to infinitely 
many energy levels 


iy+(iV, U)=pn{U), 
Z+{N,q)= ^ 


W-{N, U) = pn{U - N{N - l)/2) = p^{R), 


Z-{N,q) = 


q 


N(N-l)/2 


for 


(?)v’ ’ {q)N 

The recursive formulas obtained in the last subsection were given before in 
fermionic systems in the limit of inhnite S-values. 

If we introduce the additional energy R in place of the total energy U in fermionic 
systems the statistical weights and partition functions exhibit the same form for bosons 
and fermions. This explains the similarities previously noted, for instance the fact that 
Bose and Fermi gases (with the same number N of particles) in square-law potentials 
have the same heat capacity. Considering R instead of U enables one to let > cxo 
in the case of fermionic systems [ 0.0 0 . 0 . This limit provides with the simplest 
form of all previous formulas. 


8. Conclusion 

Non-interacting fermions and bosons statistics have been obtained in the three classical 
ensembles. The originality of the method rests on the derivation of complete and explicit 
microcanonical formulas by direct enumeration of the microstates. Canonical and grand- 
canonical relations are further derived by simple averaging. Several results are new, older 
works are referred to. 

Using the microcanonical relations requires the sole knowledge of system energy 
levels and degeneracies. Only average occupation numbers have been reported, but 
formulas for joint and higher moments may be obtained by similar arguments |l[]. 
Microcanonical formulas have been given that arise naturally. One of them leads to 
a recursive relation for the derivative of the canonical partition function. We have also 
shown that most canonical and grand-canonical quantities may be calculated from the 
one-particle canonical partition functions at temperatures T/n, n > 1. This result is 
still valid for particles that do not bear generalized charges. 

One-dimensional harmonic oscillators Schonhammer’s results p| have been 
generalized in several aspects. Bosons are also treated. Related microcanonical relations 
are given. The case of hnitely many levels is treated. Extension to arbitrary degeneracy 
shall be the object of future research. 

Other, less classical, statistical ensembles might be handled analogously. For 
example, in the grand microcanonical ensemble |^| the system may exchange particles 
with a bath. The fugacity ^ and the energy U of the system are hxed but the number n of 
its particles is allowed to fluctuate. The probability that a particular (n, f/)-microstate 
occurs is proportional to Accordingly, the grand microcanonical statistics can be 
derived from averaging the microcanonical relations in §2 with z" as a weight. 
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Appendix A. Relations between fermionic and bosonic statistical weights. 


Given the set E of the one-particle energy levels, there is a simple relation between the 
numbers W~{N,U) and W~^{N,U). It reads 

[N/2i U/2 

W-{N,U)= ^ (-l)”^hh-(n,M)fh+(A-2n,G-2M), (A.l) 

n=0 u=0 

where [iV/2j denotes the greatest integer non-exceeding A/2, and the second ^ symbol 
stands for a summation for all values of m G U between 0 and U/2. Let us sketch a proof 
of relation ( |A.1|) . Suppose the Pauli exclusion principle applies to all energy levels except 
e. Then the possible (A, [/)-con£gurations split up into those that obey strictly the Pauli 
principle and those that exhibit at least two particles at level £. Removing two particles 
at level e in the latter, denoting by W^{N,U) the number of (A, 17)-con£gurations of 
our supposed system and enumerating, result in 

W-{N, U) = hP"(A, U) - 1T"(A -2,U- 2e). (A.2) 


Further relaxing the Pauli principle at other levels one by one and iterating formula 
( |A.2|) , we get 

W-{N, U) = 5^(-l)l®liy+(A - 2\EIU - 2E,ei^^), (A.3) 


ECE 


where \E\ stands for the number of levels in E. Finally formula (|A.3|) may alternatively 
be written as formula ( |A.1|) by letting n = \E\, u = remarking that there 

are W~{n,u) ways to get u from the n elements in E. All sums terminate because 
W (A, U) = 0 for negative A or U. 

Averaging formula (|A. 1|) we readily derive the identity 


LV/2J 


Z-{N,q) = {-irZ+{N-2n,q)Z-{n,q^). 


(A.4) 


n=0 


Denoting by Z^{z,q) the grand-canonical partition function further averaging yields 
Z-{z,q) = Z^{z,q)Z-{-z‘^,q‘^), 


a fairly straightforward identity in view of (|^) . In the present case we recognize that it is 
easier to proceed the other way around, namely from grand-canonical to microcanonical 
relations. As an illustration, from the similar identities 


l = Z^{z,q)Z {-z,q), Z^{z,q) = Z {z, q)Z^{z‘^, q'^), 

we derive by expanding the series products and identifying the coefficients of z^, for 
positive A, 

N 

0 =J2i-^rZ+{n,q)Z-{N-n,q), 

n=0 

lN/2\ 

Z+(A,g)= Z-(A-2n,g)Z+(n,g2). 

n=0 
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The microcanonical counterparts of these relations are obtained by identifying the 
coefficients of 


Appendix B. Entropy 


For the reader’s convenience we link in what follows statistical quantities with 
thermodynamical ones. This can be done by examinating the expression of entropy 
in both formalisms. The von Neumann-Shannon entropy is defined as 


S = -y^p^lnp^, 

LO 


where is the probability of a state of the system and u runs through the set of all 
such states. 

In the micro-canonical ensemble with N particles and total energy 7/, the 
microstates are equally likely so that Poj(n,u) = l/kF(iV, U) and S' = In hF(iV, U), namely 
the Boltzmann entropy. 

In the canonical ensemble with N particles and temperature T, Puj(N,u) = 
q'^/Z{N,q). Thus 


S 


«eu ^ 

\nZ{N,q) + -^U{N,q). 
KbT 


qu 

Z{N,q) 


In Z{N,q) 


g In g d 
Z{N, q) dq 


Z{N,q) 


This is the classical expression of the entropy of a closed thermodynamic system with 
Helmholz free energy A{N,q) = —kBTlnZ{N,q). 

In the grand-canonical ensemble with fugacity and temperature T, the probability 
Puj{n,u) = z^q^/Z{z,q). Accordingly, we obtain as above 




u, 


Zu^u 


In 


zUqU 


MgU n>0 


Z{z,q) Z{z,q) 


= \nZ{z,q) + -^U{z,q) - -^N{z,q). 


kBT 


kBT 


The relevant thermodynamical potential (grand potential) is given by the logarithm of 
the partition function. Accordingly, entropies of subsystems add up if their partition 
functions multiply. The appropriate conditions are given in Section 6. 
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